CERN-TH/2003-070, HD-THEP-03-15 



Coherent Baryogenesis 
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We propose a new baryogenesis scenario based on coherent production and mixing of different 
fermionic species. The mechanism is operative during phase transitions, at which the fermions 
acquire masses via Yukawa couplings to scalar fields. Baryon production is efficient when the mass 
matrix is nonadiabatically varying, nonsymmetric and when it violates CP and B—L directly, or some 
other charges that are eventually converted to B — L. We first consider a toy model, which involves 
two mixing fermionic species, and then a hybrid inflationary scenario embedded in a supersymmetric 
Pati-Salam GUT. We show that, quite generically, a baryon excess in accordance with observation 
can result. 
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1. Introduction. Baryogenesis scenarios beyond the 
standard model often invoke the out-of-equilbrium decay 
of heavy particles P, 0| . Well known notable exceptions 
are the mechanisms proposed by Affleck and Dine |3| and 
by Cohen and Kaplan both being operative in the 
presence of scalar condensates. 

In this letter, we present a novel baryogenesis sce- 
nario involving nonadiabatic evolution of classical scalar 
fields during phase transitions, at which the baryon- 
lepton symmetry of the gauge group of a grand-unified 
theory (GUT) is broken. Yukawa couplings give rise to 
time dependent mass terms for matter fields and thereby 
the resonant production of particles Furthermore, 
the mass matrices are temporarily nondiagonal, induc- 
ing thus flavor oscillations between species with baryonic 
and leptonic charge. At the end of the phase transition, 
when the mass matrix becomes diagonal, the charges get 
frozen-in. Since this mechanism relies on the interplay of 
coherent particle production and (_B— L)-violating flavor 
oscillations we call it coherent baryogenesis. We empha- 
sise the conceptual simplicity of coherent baryogenesis, 
since it involves the tree-level dynamics of quantum fields 
only. 

2. Formalism. For our calculations we adapt the ki- 
netic approach (c/. Refs. which makes use of the 
dynamics of Wigner functions, which arc closely related 
to classical phase space distributions. 

To compute them, wc extend the formalism of Ref. 7, 
0. We define the Wigner function. 



d'^re'''-'' {Q\Mx - rl2)i,a{x + r/2)|0). 



where a,b — 1, ...,N are species indices and (^7*^5^)^ = 
ij^S^ is hermitian. Our derivation will go through for a 
2-point function with general density matrix, but in view 
of our applications in inflation we prefer to write it with 
respect to the vacuum from the outset. Introducing an 



N X N matrix M and its hermitian and antihermitian 
parts, respectively, 

Mff = i (Af + Mt) , Ma = y^ {M - Mt) , 
we find that iS^ obeys the Wigner space Dirac equation 



+ ^7°9t - {Mh + il''MA)t 



\dtdka 



(1) 



The mass matrix M emerges generically from Yukawa 
couplings to scalar field condensates, £yu = —y4>4'R'4^L + 
h.c, which can induce CP-violation in the fermionic 
equations 0, O, 0, 0|. It is noteable, that in coher- 
ent baryogenesis the condensate does not have to carry a 
charge, while this is necessary for the Affleck-Dine mech- 
anism ^3], therefore being conceptually different. This 
can be seen explicitly from the two examples in this let- 
ter, where the scalar condensates (j) are always real, and 

therefore = dt (j)) = 0. 

However, the time dependence of M plays here another 
important role. For = 1, the Dirac equation consists 
of two first order differential equations due to the two de- 
grees of freedom of the fermionic field. Thus, also dM/dt 
can contribute CP-violation, and in general, both sources 
from AI and dM/dt cannot be simultaneously removed 
by local phase reparametrizations of the fermionic fields, 
cf. Ref. ^8] for an example. When A^ > 1, even higher 
derivatives of M are involved, allowing, in principle, mul- 
tiple sources of CP- violation. We stress that this is a very 
different situation from the Standard Model quark mix- 
ing Cabibbo-Kobayashi-Maskawa matrix, or lepton mix- 
ing Maki-Nakagawa-Sakata matrix, where at least three 
generations of quarks or leptons are required for one CP- 
violating phase. 

We now make use of the fact that the helicity operator 
/i = k • ^O'y^'^ , commutes with the Dirac operator in 
and decompose the Wigner function 
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ijos;: = -{l + hk- 



(2) 



where we have omitted the species indices, k — k/|k| 
and cr'', (/i = 0,1,2,3) are the PauH matrices. We 
multiply by p'^, take the trace and integrate the her- 
mitian part over feg- Introducing the 0th momenta of g^^, 
ffih — J idko/2Tr)g^h, we obtain the system of equations 

foh + ^[MHJlh]+^[MAJ2h] = (3) 
fih + 2h\k\f2h + i [Mh, foh] - {Ma, hh} = 
f2h-2h\k\hh + {M„,hh} + t[MAjoh] = 

f3h~{MH,f2h} + {MAjlh} = 0. 

The /^ft(a;,k) can be interpreted as follows: foh is the 
charge density, /a/,, is the axial charge density, and fih 
and f2h correspond to the scalar and pseudoscalar den- 
sity, respectively. Note that the commutators in JSJ, 
which mix particle flavors, are essential for the produc- 
tion of the charges foh, and thus for our scenario. 

For an originally diagonal mass matrix and adiabatic 
conditions, the initial Wigner functions describing zero 
particles are (c/. Ref. Ql): 
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with 
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where LOa = \/k^ + |Afaap. For a nondiagonal, but her- 
mitian, M , one obtains initial conditions by an appro- 
priate unitary transformation. If additionally Ma 0, a 
biunitary transformation is necessary for diagonalisation. 

Since /oh is the zeroth component of the vector cur- 
rent, the charge of the species a carried by the mode 
with momentum k and helicity h is simply qah{k) = f^^. 
From cqs. Q one can derive that, in order to generate a 
nonvanishing charge qah{k), i.e. f^^ 0, M should not 
be symmetric. Note also, that the Lagrangean 

C = i^a^a - MMh + ij'^MA)bai^a 

is U{1) symmetric, and thus Qahik) is conserved. 

3. Toy model. We now consider a two species model, 
where fundamental SU (2) fermions couple to an adjoint 
scalar triplet and to a singlet such that the mass 
matrix is given by M = 4>°1 -I- $V. While we fix the 
fields associated with the diagonal generator as = /i 
and = ijl/2, with /j, being a mass scale, we let and $^ 
move freely in a harmonic potential, starting form arbi- 
trary initial conditions. Conjugating the fermionic sector 
of ||2Jl under CP while leaving the scalar condensate in- 
variant, we find that, in spite of M being hermitian, CP 



is broken for the fermions, because M ^ M* . Damp- 
ing is introduced through a phenomenological decay rate 
F and through the Hubble expansion in a matter domi- 
nated universe, e.g. the scale factor is a = OmTf' , where 77 
denotes the conformal time, and ' = d/drj — ad/dt. The 
equation of motion for a scalar $(77) is given by 



„ a' I n dV I 
$" + 2—$' + a^—+ aF$' = 0. 



a 

Writing bj\^d'^V / d^'^ and setting F = 0, the solutions are 



(5) 



Cl 



■ cos 



(6) 



For we employ ci 
Cl = 0, C2 = /i and w$ 



/I, C2 = and uj^ = /x, for $ , 
1.5/i, and we set am = is^. 
We approximate the effect of damping by multiply- 
ing the solutions ^ by A = exp (— ^OmFTy'^) , where 
F — 0.1/x <C The equations of motions for the 

Wigner functions in conformal space-time are then sim- 
ply obtained by replacing M by aM in Q [3, H- We 
illustrate the motion of the mixing fermionic mass terms 
in conformal time in figure ^ 
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FIG. 1: Parametric plot of the motion of a{rj)A^i and 
a{Tj)A$2 for ri G [2.3^"', 4^1-^]. 

Requiring that there are no particles at 77 = 2.1/i^^, 
we choose initial conditions in accordance with 10} and 
solve (jSl numerically, to find fermion number production 
as displayed in figure El as a function of the conformal 
momentum k. We sum over the helicities, qa = qa+ + qa~, 
and note, that in the present case qa+ = qa-, for Ma = 0. 
Integration over fc, qa = f dkk^qa{k)/2TT^ , gives the 
charge densities qi — —q2 = 2.4 x 10~^(/i/a)^. If qi was 
charged under B, our toy model would lead to successful 
baryogenesis. 

4. Hybrid inflation in a supersymmetric Pati-Salam 
model. We shall now discuss the implementation of co- 
herent baryogenesis in a more realistic model. In or- 
der to generate a baryon asymmetry which survives the 
sphaleron washout, we require the presence of {B — L)- 
violation. This is the case in several GUTs, e.g. E{6), 
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FIG. 2; The fermion numbers qi+{k) and q2+(k) in the toy 
model. 

a subgroup of which is S'O(IO), and in the Pati-Salam 
group, Gps = SU{'i)c X SU{2)l x SU{2)b., which ap- 
pears as an intermediate stage of breaking of 5*0(10) 
down to the Standard Model. For simplicity, here we 
study an extension of the hybrid inflationary scenario 
embedded in a supersymmetric Pati-Salam model, which 
is considered in jOj and does not suffer from the monopole 
problem. The relevant terms of the superpotential are 

2 



Ms 



(7) 



Under Gps, the fields transform as H'^ = (4, 1, 2), H'^ = 
(4,1,2), S= (1,1,1), G= (6,1,1) lOj. We adopt the 
notation 



''HI 



'■H2 



JC JC JC 

"-HI "'H2 "H3 



H3 '^H 



and likewise for . With SU{A)c broken to the Stan- 
dard Model, G = D+D, with D = 3 and D = 3 oi SU{3)c. 
Vanishing of the D-terms requires H'^ = H" , which we 
assume throughout this discussion. During inflation, the 
neutrino-like component vfj of H"^ has the vacuum ex- 
pectation value 



{Wj,\) = Ms[n/i2P)Y- 



(8) 



which evolves during the waterfall regime to the super- 
symmetric vacuum, where S = and 



{Wh\) = 



2(3 



(9) 



All other components of H'^ as well as the fields G vanish 



at all times. Thus, IH"]^ = H'^H" = 
implies the following scalar potential: 



V=2 
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, and eq. 
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Numerically, we simulate the evolution of the scalar 
fields from the end of inflation until they settle to the 



K = 0.007 




= 1 


c = 


0.01 


F = 0.008^ 


Ms 




c = 


0.01 xexp (ixlO^^) 



TABLE L Parameters used in numerical simulation. 

supersymmetric vacuum. The choice of our parameters 
is listed in table for a discussion of their implications 
for cosmic observations, see Ref. 0|. We let point 
in the real direction and remark, that not all constants 
can be made real by field redefinitions, thereby allowing 
a nontrivial phase between C and ^. Note, that the sub- 
stitution C < — > ^* results in the opposite sign oi B — L 
and B and hence, when C — there is no generation of 
B—L. Furthermore, that Ms is rather small is not an es- 
sential feature of the model, but was introduced to avoid 
large initial Higgs field amplitudes, which are unplesant 
to deal with numerically. 
3 




FIG. 3: Evolution of the field ;/|r- 

Eventually, condenses to the value given by (c/. 
figure|3I), and the field S settles to zero. In this situation, 
the flavor oscillations between the Dirac fermions 



and X2j = 7 

are scotched, since their mass term, 

( m^H)i] irn^ \ 



(10) 



( Xlj X23 ) 



where md—{S) (k/2 — /3(|t/|f |)^/M|) , becomes diagonal. 

Starting with initially zero fermions, we compute the 
charges qij{k) and q2j{k) associated with xij s-iid X2j, 
respectively, at the time when the phase transition is 
completed, by using ©, as shown in figure^] Integration 
then yields qi — —92 — 1-0 x 10~^(^/a)'^, where we have 
already summed over the three colors, and where in our 
computation a w 120 during the phase transition. 

In the presence of the {B — L)-violating condensate of 
the right handed neutrino i^'^ acquires a Majorana 
mass from the additional, nonrenormalizable contribu- 
tion jP'^H'^F'^H'^ /Ms to the superpotential {Tj), where 
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FIG. 4: The produced charges of the Dirac fermions xij, X2j, 
summed over both hehcities. 

F'^H'^ form a gauge singlet, and F"^ — (4,1,2) are the 
superfields containing the right handed quarks and Icp- 
tons, 

pc ^ ( "r ul \ 

\ di di e= y ■ 

This also gives rise to the Lagrangean terms 

l{{l^H)lMs)[^V=(l)d':'lpd''^. + i^d-fj^u-fpd'l] +C.C., 

allowing the decay of the c?^^.-component of xij in H10|l 
to d''*+ v'^* , where one of the latter particles is fermionic, 
the other scalar. The Majorana neutrinos v'^ are their 
own antiparticles; neglecting the small effects induced by 
possible mixing angles and a CP- violating phase in the 
Majorana mass matrix familiar from leptogenesis scenar- 
ios |2l ll4lll6lll7l |. their decay leaves behind no net charge. 
Note that in turn, our mechanism does not require any 
lower constraint on these angles and the phase as in lep- 
togenesis. 

The coupling H" F" H"" / M s allows the decay of the 
ff'^-fields through the reaction from the c?^-componcnt 
of X23 in l(Tn|l to d'^+ u'^. The charges hence get trans- 
formed to (_B — L) = — |g2- This number is promoted 
by sphaleron processes to _B = (10/31)(i?— L) loj, where 
we assumed two complex Higgs doublets. 

Hence, the final value oi B — L arises here due to the 
transformation of other charges in decay processes. How- 
ever, we point out, that models are conceivable where 
coherent particle production directly leads to standard 
model particles, the (i? — L)-charge of which is conserved 
in the subsequent history of the universe. 

To obtain a lower bound on the baryon-to-photon ra- 
tio, we need an upper bound for the entropy, which can 
be obtained by assuming a complete and instantaneous 
conversion of the vacuum energy g— [K^Af|/(4/3) — Kji^Y 
into a thermal bath of highly relativistic particles, with 
the energy density g = 7r^g*r**/30, where g* = 221.5 is 
the number of degrees of freedom in the MSSM. The 
entropy per unit volume is then s = 27r^g*T'^/45, and 



we thus estimate the generated baryon-to-photon ra- 
tio to be B /n^ ~ 9.4 x 10"^". Also, a contribution 
from boson production may arise, which we do not con- 
sider here. Note that our result was obtained by choos- 
ing a small CP-violating phase (c/. tabled, indicating 
that there is an ample phase space of couplings, which 
leads to baryon production consistent with the observed 
B/n^ = 6.1 ± 0.3 X 10-1° In conclusion, we have 

demonstrated that coherent baryogenesis is a viable, effi- 
cient and natural candidate for the creation of the baryon 
asymmetry of the universe. 
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